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We review two definitions of temperature in statistical mechanics, Tb and To, corresponding to 
two possible definitions of entropy, Sb and So, known as surface and volume entropy respectively. 
We restrict our attention to a class of systems with bounded energy and such that the second 
derivative of Sb with respect to energy is always negative: the second request is quite natural 
and holds in systems of obvious relevance, i.e. with a number N of degrees of freedom sufficiently 
large (examples are shown where N 100 is sufficient) and without long-range interactions. We 
first discuss the basic role of Tb, even when negative, as the parameter describing fluctuations of 
observables in a sub-system. Then, we focus on how Tb can be measured dynamically, i.e. averaging 
over a single long experimental trajectory. On the contrary, the same approach cannot be used in 
a generic system for To, since the equipartition theorem may be spoiled by boundary effects due to 
the limited energy. These general results are substantiated by the numerical study of a Hamiltonian 
model of interacting rotators with bounded kinetic energy. The numerical results confirm that the 
kind of configurational order realized in the regions at small Sb, or equivalently at small |Ts|, 
depends on the sign of Tb. 

PACS numbers: 


I. INTRODUCTION 

Two different definitions of temperature in equilibrium statistical mechanics have been recently the subject of an 
intense debate [Ml , after the publication of experimental measurements of a negative absolute temperature PH. 
In m it was demonstrated the possibility to prepare a state where the observed distribution of the modified kinetic 
energy per atom appeared to be inverted, i.e. with the largest population in the high energy states, yielding a de facto 
negative absolute temperature. 

The possibility of a negative absolute temperature is well known since the theoretical work by Onsager on the 
statistical hydrodynamics of point vortices [l3| and the experimental and theoretical results on nuclear spin systems 
by Pound, Ramsey and Purcell (see [l^ - [l^ for a review and discussion). In those investigations, it was clear that an 
inverse temperature parameter /3 ranging in the full infinite real line (—oo, oo) did not lead to any inconsistency or 
paradox. Ramsey in 1956 already realised that “the Caratheodory form of the second law is unaltered.” 0 

A negative absolute temperature appears whenever the microcanonical entropy is non-monotonic in the energy, a 
condition which can be realized when the total energy has a global maximum, which may happen when the phase space 
is bounded. There are also cases where the phase space is bounded but the energy diverges: again this may lead to a 
non-monotonic entropy, an important example is given by point vortices [131 I17 h 21| . It is crucial to highlight that the 
lack of monotonicity (for entropy vs. energy) is realised if one adopts the simplest definition of microcanonical entropy, 
which is related to the logarithm of the number of states with a given energy. Since such a definition appears in the 
so-called “tombstone formula” written on Boltzmann’s grave, “S = fclogIT”, it is often referred to as Boltzmann’s 
definition of entropy. Even if not historically precise [8|, we adopt the same convention (but setting k = 1) and 
call “Boltzmann entropy” of a system with Hamiltonian H{Q, P) - where Q and P are vectors in , being d the 
dimension of the system - 

SB[E,N)=\ogu:{E), (1) 

being uj{E) the density of states, i.e. 

w{E) = JS{H- E)d^^Qd‘^^P = (2) 
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and S(£') the total “number” of states with energy less or equal then E, that is 


^{E) = f ( 3 ) 

J H<E 

In definition o we have ignored an additive constant which is not relevant in our discussion. In Q it is stated that 
the validity of the second principle of thermodynamics depends on the value of this arbitrary constant. Nonetheless, 
such an arbitrariness and the consequent paradox can be removed if all the quantities (energies, positions, momenta, 
time etc...) are considered adimensional. Propagating the denomination, it is customary to define the “Boltzmann 
temperature” through 



dSeiE, N) 
dE 


( 4 ) 


Some authors [I|, Q have argued that a different definition of microcanonical entropy, proposed by Gibbs, has to 
be used in statistical mechanics, in order to be consistent with a series of “thermodynamic” requirements and avoid 
unpleasant paradoxes. The Gibbs entropy, which is always monotonically increasing, reads 


SG{E,N)^\ogi:{E), 


and leads to the Gibbs temperature definition, which is always positive: 


Pg = 


1 


dScjE^N) 

dE 


( 5 ) 

( 6 ) 


Let us note that, since Tb is defined directly on the surface of interest (i.e. that at constant energy E), from the 
point of view of the ergodic approach its use appears rather natural. The Gibbs temperature, on the other side, enters 
through an ensemble average in the equipartition formula of textbooks : 


) ( 7 ) 

where Xi is any of the components of vector (Q,P) and the average is done in the microcanonical ensemble. In 
Section m we will discuss the limits of application of formula 0 when the energy is bounded. We also mention 
that Tg appears in the theory of Helmholtz monocycles (which had an important role in the development of the 
Boltzmann’s ideas for the ergodic theory), for one-dimensional systems 1^. 

In spite of the fact that, in our opinion, the basic features of the different definitions of temperature do not present 
particular technical or conceptual subtleties, there is a certain confusion in the literature; therefore a general discussion 
of the topic can be useful. In this paper we present a line of reasoning where Boltzmann temperature Tb (positive or 
negative) is the (unique) proper parameter which is relevant for the statistical properties of the energy fluctuations, as 
well as in determining the flux of energy between two systems at different temperatures, in addition it is measurable, 
without the appearance of any evident inconsistency. Let us remark that the systems discussed in Q , from which the 
authors try to show that only Tg is the “good” temperature, are small {N = 0{1)) and/or with long interactions. 

In Sectionini after presenting the class of physically relevant systems which are the subject of our study, we describe 
how the Boltzmann temperature Tb naturally describes fluctuations of observables in subsystems, in analogy with the 
derivation of the canonical ensemble from the microcanonical one. In Section mi we discuss dynamical (“ergodic”) 
measurements, which can reproduce Tb but are in general unsuited to measure Tg' in particular we show a possible 
failure of the equipartition theorem. In Section llVI we report a series of numerical results with a model of interacting 
rotators with bounded kinetic energy, discussing the many practical uses of Boltzmann temperature. Summary and 
conclusions are drawn in Section El together with a critique of some of the arguments used, in Q, to rule out the 
thermodynamic meaning of Tb- 


II. THE RELEVANCE OF THE BOLTZMANN TEMPERATURE 

In this section we show, following the standard approach that can be found even in some textbooks, the unavoidable 
role of Tb in many problems of statistical mechanics. 
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A. Systems of physical relevance 

In the rest of the paper we consider systems made of a finite but large number iV 3> 1 of particles with local 
interactions, i.e. we exclude long-range potentials or mean-field models. It should be understood that long-range 
interactions certainly widen the phenomenology of statistical mechanics and may lead to complicate functional de¬ 
pendences for Sb{E, N), e.g. with several maxima or minima, even for large N. Nevertheless they are not necessary 
for the discussion of negative temperature and, most importantly, they represent quite a peculiar case where even 
thermodynamics is not obvious: for instance, it is not evident that the typical Gedankenexperiment of putting in 
contact two - previously isolated - systems can be realized, as the isolation condition is prevented by the long-range 
interaction. 

We also assume that Sb{E,N) is always convex, i.e. d'^SB{E,N)/dE^ < 0. This is certainly true in the limit 
of vanishing interaction and in short-range-interacting systems for large N, since S's is strictly related to the large 
deviation function associated to the density of states . Let us stress that these large values of N are not necessarily 
“thermodynamic” {N —)• oo): for instance in Sec. IIVI we will exhibit a system that possesses all the required features 
already at N = 100. In general such a value of N will depend on the specific system, corresponding to situations 
in which some common approximations (e.g. Laplace approximation for exponential integrals) can be safely applied. 
In Sec. 11.C we discuss in some details the origin of the convexity of Sb{E,N). It is easy to understand that this 
assumption implies the validity of the second principle of thermodynamics, as discussed in the next subsection. 


B. Second law and energy flux between two systems in contact 

Let us consider a system A of particles described by the variables {Q^,P^} and Hamiltonian LIyt(Q_ 4 , P^), 
a system B of Nb particles described by the variables {QB,Pe} and Hamiltonian LIb(Qb,Pb) and a small coupling 
among the two, so that the global Hamiltonian is 

H = HAQa, Pa) + Hb{Qb, Pb) + Qb)- (8) 

If the two Hamiltonians have the same functional dependencies on the canonical variables (i.e. they correspond to 
systems with same microscopic dynamics, with possibly different sizes and Nb), for large N, we can introduce 
the (Boltzmann) entropy per particle 


SB{E,N) = NS{e) 


E 

iV’ 


(9) 


with S{e) a convex function, identical for systems A and B. Let us now suppose that systems A and B have, 
respectively, energy Ej^ = and Eb = Nb^B and the corresponding inverse Boltzmann temperatures and 

When the two systems are put in contact, a new system is realized with N = N_a + Nb particles. Let us call 
a = Nj^^/N the fraction of particles from the system A. We have that the final energy is Ef = Ej^ + Eb = Ncf, where 
e/ = ae_A + (1 — a)eB and final entropy 


SB{Ef,N) = NS{ef) > NAS{ei) + NeSiee) = N[aS{eA) + (I - a)S{eB)]. 
The previous inequality follows from the convexity assumption for S{e) which implies 


( 10 ) 


S{aeA + (1 - a)eB) > aS{eA) + (1 - a)5'(eB). (11) 

The final inverse temperature is intermediate between and e.g. if ea > - that is - then 

. ( 12 ) 

The energy flux obviously goes from smaller Pb (hotter) to larger Pb (colder). The consequence of convexity is that 
Pb{E) is always decreasing and a negative value does not lead to any ambiguity. Confusion may arise from the fact 
that Tb < 0 is, for the purpose of establishing the energy flux, hotter than Tb > 0. However if Pb is used, the 
confusion is totally removed M- 

We also briefly discuss a particularly interesting case with different Hamiltonians. Suppose that for the system A 
negative temperatures can be present, whereas system B has only positive temperatures; it is quite easy to see that 
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the coupling of the system A at negative temperature with the system B at positive temperature always produces a 
system with final positive temperature. Indeed, at the initial time the total entropy is 

Si = S^{Ea) + S‘^{Eb), (13) 

while, after the coupling, it will be 

SF = S^{EU) + S’^iE'^), (14) 

where E'^ + E'j^ = E^ + Es and, within our assumptions, E'^ is determined by the equilibrium condition that 
Sf takes the maximum possible value, i.e. 


dE'^ dE'j^ ■ 


(15) 


Since /3g is positive for every value of iSlg, the final common temperature must also be positive. The above conclusion 
can also be found, without a detailed reasoning, in some textbooks [^.1^. 


C. Subsystems 


Let us consider a vector X in (with < N), that is a subsystem of the full phase space (Q,P), and let us 

indicate with X in remaining variables. We have 


H = iJi(X) + iL2(X) + Hi{X, X) 


(16) 


with an obvious meaning of symbols. 

Let us consider the case N ^ 1 and Ni N. In the microcanonical ensemble with energy E, the probability 
density function (pdf) for the full phase space (Q, P) is 


P(Q,P) = 


1 


j{E,N) 


S{H{Cl,P)-E). 


(17) 


The pdf of X can be obtained from the latter, by integrating over X. If the Hamiltonian iL/(X,X) is negligible (a 
consequence of our assumption for non long-range interaction) then we have 


P(X)c. 


uj{E-Hi{X),N-Ni) 


u}{E,N) 

It is now possible to exploit the definition of Sb and get 

which, together with (fTS|) leads to 

P(X) oc 


(18) 

(19) 

( 20 ) 

( 21 ) 


When Ell is bounded (as in our assumptions), the previous simple derivation can be done irrespective of the sign of 
Pb- It is immediately clear from the above argument that Tb is the temperature ruling the statistics of fluctuations 
of physical observables in a subsystem. For instance, the pdf of the subsystem (i.e. the canonical ensemble) energy 
El reads 


P(F;i,iVi) oc uj{Ei,Ni)e-^^^^ ^ ^ISb(EfNi)-i3bE,]_ 


( 22 ) 


Of course the above result holds in the (important) case where the two subsystems are weakly interacting and Hi E. 
Therefore, for ei = Ei/Ni, one has 


P(ei, A^i) oc 


(23) 


which is a large deviation law where the Cramer’s function C{ei) is C(ei) = Pb^i — S{ei) + const. From general 
arguments of theory of probability, we know that - if a large deviation principle holds - > 0 so ^ < 0. 

The validity of the large deviation principle can be easily shown for non-interacting systems. For weakly interacting 
systems it is quite common and reasonable, and can be stated under rigorous hypothesis 
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D. The generalised Maxwell-Boltzmann distribution 

The extreme case of the above considerations is when = 1, that is to say the fluctuations of a single degree 
of freedom (e.g. a momentum component of a single particle) are observed. This becomes interesting when the 
Hamiltonian has the form 


N N 

H = ^g{Pn)+^V{qn,qk) (24) 

n—1 n,k 

where the variables {pn} are limited and the same happens for the function g{p). 

Repeating the arguments in the previous subsection, one may compute the probability density for the distribution 
of a single momentum p, obtaining 


P{p) 


^{E-g{p),N-l) 

uj{E,N) 


(25) 


which, again, is valid for both positive and negative 13b- We mention that in the experiment in [TTI |. the above recipe 
has been applied to measure both positive and negative system’s temperatures. 

From Eqs. (1221) and (1251) the true deep meaning of the (Boltzmann) temperature is quite transparent: it is a quantity 
which rules the pdf of energy of a subsystem (or the momentum of a single particle). Let us note that since Tb is 
associated to the large microcanonlcal system (In physical terms the reservoir) it is a non-fluctuating quantity [29| 
also for each sub-system and, in general, for non-isolated systems. In the conclusions, we discuss again such an aspect 
which is not always fully understood, see e.g. Ref. Q 


E. Temperature and order 

In usual statistical mechanics, low temperatures - or, better, high values of Inverse temperature - are usually 
associated to the possibility of some kind of order, the most noticeable example given by phase transitions. Intuitively, 
one would expect such a situation whenever uj{E) Is relatively small, which usually corresponds to regions where |/ 3 b | 
is large irrespective of the temperature’s sign. A famous example where such an order at negative (small) temperatures 
was observed Is that of polntllke vortices discussed by Onsager in [T^ . The system, obtained as a particular limit 
from two-dimensional Euler equations, describes N points of vorticities {Fi, ...,F 7 v} in a two-dimensional domain fl: 
the equation of motions of the coordinates (Xn, Un) of the n-th point vortex are shown to be (see for instance 0) 

" dt dy^ ' " dt dx^ ^ ^ 


with Hamiltonian 


H = (27) 

where t/(r) is the Green function of the Laplacian in H: in the infinite plane one has G{r) = —l/47rlnr where 
Xij = ^{xi — Xj)'^ E {yi — ViY- The canonical variables In this case are 

qi = > Pi = •\/il^sign(Fi) yi (28) 

Onsager showed that If the domain of H is bounded, then negative Tb are achieved at large values of the energy. At 
large energies a particular spatial order appears too: clusters of vortices with the same sign of the vorticity are the 
structures most easily found. It is interesting to notice that Tb < 0 (and the corresponding clusterization) is not 
a peculiarity of the divergence of Q{r) in r = 0, nor of the long range nature of the interaction: indeed, it can be 
obtained with any arbitrary t/(r) having a maximum (even finite) in r = 0, and vanishing at large r, provided that 
the domain is bounded. The presence of spatial order at high values of energy, in the form of discrete breathers, has 
been observed also In the discrete non-linear Schrddinger equation and analogous systems [13, [m. In Section 113 we 
introduce a different, in a way simpler, model which still exhibits spatial order at small negative temperatures. 
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III. HOW TO MEASURE Tb AND Ta 

The definitions of /3 b and Pc given in Eqs. o and m are based on the functional dependence of the phase space 
occupations oj{E) and S(i3) upon the energy. In a real or numerical experiment it may be cumbersome or even 
impossible to make use of those definitions to measure the two temperatures: for instance, an empirical estimate of 
oj{E) (and therefore of S(i3)) will always be limited by the available statistics (number of independent measurements 
of E) and therefore cannot provide a clear answer, for both /3 b and Pc^ in ths interesting regimes where u}{E) ^ 0. 

On the other hand it has been shown that /3b can be obtained as a microcanonical average of a certain 
observable. The recipe is the following 

/3b =< i?(X) > , i?(X) = V • (29) 

where V stands for the vector of derivative operators along the degrees of freedom in the full phase space X = (Q, P). 
From (0^ one has, assuming the ergodicity, that /3b can be computed with a molecular dynamics simulation, and, at 
least in principle, by a long-time series from an experiment. It is interesting to notice that such a kind of recipe does 
not exist for Sb{E,N) or Sa{E,N) [s^- It is clear that, in view of the considerations in Sections III Cl and III D1 one 
may always measure fluctuations of appropriate observables, such as subsystem’s energy or single particle momentum, 
to get an estimate of Tb. 

Coming to /3g, a way, even discussed in textbooks and considered sometimes rather important to approach the 
problem of its measurement is via the equipartition theorem, which states 

However the usual derivation of Eq. (15(11) implies the possibility to neglect boundary terms in an integration by parts. 
Such a possibility is challenged in the class of systems with bounded energy and phase space that we are considering. 
In particular it is easy to show that (1301) does not hold under the simultaneous realization of the following conditions: 

• bounded space of the canonical variables; 

• bounded derivatives of the Hamiltonian 

• bounded energy from above and below: E^ < E < Em] 

• vanishing density of states at the boundaries, i.e. ui{Em) = 0. 

Given such conditions, one has that, on one side. 


Tg{E) 


E{E) 

uj{E) 


(31) 


diverges when E —)• Em- On the other side, is limited, resulting in a contradiction. 

A failure or the equipartition formula Eq. (1301) is also possible in systems where there are no negative temperatures, 
i.e. Tq ~ Tb > 0 for all E. Consider, for instance, the following Hamiltonian 


H ^ Y-f e^(l - cos((()„ - 0„_i)) (32) 

n—1 n—1 


where G [—tt, tt). For large T, i.e. E ^ eN, the contribution to T{E) of the variables {pn} does not depend too 
much on the value of T, so that 

T,{E) ~ Eo(T) oc E^^^ , (33) 

and Tq ~ 2E/N and, for large N, Tb = Tq + 0{1/N). 

On the other hand it is easy to see that 


dH 

'‘d(j)n 


< 2'Ke 


and, therefore, the equipartition formula {(f)n§^) = Tq does not hold for large value of E and N. 


( 34 ) 
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FIG. 1: Phase space sampling: we report the reconstruction of the density of states uj{E) and its integral E(i5) = dE'uj{E'). 
The two functions are normalized with 'E{Em = 2N{1 + e)). The parameters of the system are: N = 100 and e = 0.5. 


IV. NUMERICAL RESULTS FOR A SYSTEM WITH NEGATIVE TEMPERATURE 

In this Section we present a detailed study of a system composed of N “rotators” with canonical variables 
(j)i, ,Pi---Pn with all (j)i and pi defined in [—7r,7r), and with Hamiltonian 

N N 

H{4>1,.. .,(j)N,Pi, ■ ■. ,pn) = “ cos(p„)] + e^[l - cos(0„ - 0„-i)]. (35) 

n—1 n—1 

Choosing, as boundary condition, t/fo = 0 guarantees that the only conserved quantity by the dynamics is the total 
energy E. The equations of motion for the rotators can be readily obtained applying Hamilton’s equations to Eq. 

( 1551 ) : 


= sin(p„), 

Pn = -e + sin((()„ - (j)n+i)) ■ (36) 

It is immediate to verify that the energy has a maximum value Em = 2V(1 + e) which is realised when Pn = and 
(/)n - 4>n-i = TT for every n. 

When e = 0 it is immediate to see that Hamiltonian in Eq. (1351) implies negative Boltzmann temperatures. Indeed 
at small energy one has 1 — cos(p„) ~ p^/2 so that 


T.{E) - CnE^/^ , uj{E) ~ (37) 

with Cjv = r(N/ 2 +i) ■ Close to Em = 2N one has 1 — cos(p„) ~ (tt — pn)^ 12, therefore when E approaches Em 

it is 


N 


EiE) = E(Em)-(2tt)^ f l[dpnc^ E(£;M)-(27r) 

JE<H<Em n—1 


N 


N 




, n dpn = ^(Em)-Cn(Em-E)^/^ 


( 38 ) 


and therefore 


u:(E) ~ ^Cn(Em - E)^/^-^ . (39) 

In conclusion we have that (jj(E) = 0 if E = 0 and E = which implies a maximum in between and a region 
(at high energies) with negative Pb- The previous scenario is expected to hold also in the prescence of a small 
interaction among the rotators and can be numerically confirmed with a sampling of the phase-space (see Fig. [IJ: 









random configurations of the system are extracted with an uniform distribution over the phase space and is 

reconstructed by counting the number of configurations lying in a small interval of width 5E around the energy E. 
It is rather evident from Fig. [T]that: the density of states u^E) has a maximum in i? « EmI^\ it is an increasing 
function for E < E whence Tb > 0; it decreases ior E > E whence Tb < 0. Unfortunately, such a sampling is reliable 
only in a narrow region around E: indeed, there are very few configurations with energies much larger or smaller than 
E and, therfore, there is an extremely small probability to extract such configurations with this procedure. 

For this reason, we have performed dynamical measures through numerical simulations of the motion of the system: 
the integration of Eqs. (I5S1) is done with the usual Verlet scheme with a time step At = 10“^. 


A. Measure of Tb 


Measurements of the Boltzmann temperature are done with the two methods discussed in the previous Sections. 
In particular, by computing the following average (over a single trajectory of the system) 


N 


p{p) = lim ^ / dt 5 {pi{t) - p), 
r^-oo At 7o ^ 


(40) 


for different values of p, and assuming that the system is ergodic, we recover the single-particle-momentum probability 
density function P(p), Eq. (1^ . The result of such a measure is reported in Fig. [5] for two different values of 
energy E+ < E and E- > E the measured p{p) is plotted as a function of the “kinetic energy” of the individual 
rotator g(j)) = 1 — cos(p). The presence of a negative temperature at E = can be readily indentified by means 
of the consideration in Section III D1 Indeed, on one hand, the exponential behaviour of p(p) guarantees that the 
approximation used to obtain Eq. (1^ is already valid (for every value of g{p)) at N = 100. On the other hand, the 
clear positive slope of the function at E = is a direct consequence of the fact that Tb{E-) < 0: the opposite 
situation is encountered at E = E+, where the decreasing behavior of p(p) indicates a temperature Tb(E+) > 0. 
These conclusions can also be drawn by measuring the time average of the function R{X), Eq. (I^Hl) : in the inset of 
Fig. [^we report the temperature obtained with the cumulated average of i?(X) up to time t, namely 

5^ = W)). (41) 

for E = i?_|_ and E = E-. These two quantities converge, for large t, to an asymptotic value representing an estimate 
of the inverse Boltzmann temperature /3b of the system. This value, as expected, is positive for E = E^ and negative 
for E = E-\ moreover, the values are in very good agreement with the slopes of the single particle distribution 
function, as shown by the dashed and solid lines in Fig. [51 


B. Equivalence of ensembles and the equipartition formula 

Let us briefly discuss the problem of the equivalence of ensembles. In the usual treatment of textbooks one starts 
from Eq. (23): assuming that S(e) is convex and performing a steepest descent analysis, for large N, one obtains the 
canonical functions from the (Boltzmann) microcanonical ones, e.g.: 

TB{e)S{e)=e- f{TB{e)), (42) 

where f{T) is the free energy per particle in the canonical ensemble. In addition the energy fluctuations are negligible. 
In such a derivation, the relevant point is only the convexity of S{e) and nothing about its first derivative is asked. 
Therefore, the equivalence of ensembles naturally holds under our hypothesis even for negative Tb- Since Tb and Tq 
can be different even for large iV, as in our model defined with Eq. (35), it is evident that Tg is not relevant for the 
ensemble equivalence. 

A common way @ to measure the Gibbs temperature is by means of the equipartition formula, Eq. (1301) : for the 
Hamiltonian in Eq. dSSl) one should get 


{pk s\npk)E = Tg{E), (43) 

for every 1 < k < N. In the present subsection, we use the notation ()e to denote the average in the microcanonical 
ensemble, in order to distinguish it from a canonical average Op which is useful to get some analytic expressions and 
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FIG. 2: Measure of the Boltzmann temperature in the rotators chain for N = 100 and e = 0.5. Probability distribution 
function of the momentum of the rotators as a function of their “kinetic energy” g{p) = 1 — cos(p) at energy E = E- — 170 
(blu squares) and E = E+ = 130. The slopes of the two black straight lines are l/T"^ [E), where T^{E) is the asymptotic 
value of the corresponding curve in the inset. Inset: The Tr obtained from the cumulated average of the observable i?(X(f)) 
over a trajectory up to time t a-t E — 170 (blue line) and E = 130 (red line). 


better investigate the validity of Eq. (P)) . The canonical probability density reads 

p(0i,... ... ,p„) = (44) 

where Z{f3) is the partition funcion and /3 the (external) inverse temperature, that can be either positive or negative: 
if such a distribution is derived from a larger isolated system, as already discussed in Section Hi Cl the temperature 
in the canonical ensemble is precisely the Boltzmann temperature of the whole system. A simple explicit expression 
(see details of analogous calculations in Ref. [s^) can be derived for the mean energy 




^hWe) \ 

We) j ’ 


(45) 


where Io{x) and Ii{x) are, respectively, the zeroth and the first modified Bessel function of the first kind. Analogously, 
one can get an analytic formula for the equipartition function 


{psin{p))p 


1 e ^ 


(46) 


Let us remark that Eqs. (HSI) and (H51) hold for both positive and negative /3. In Fig. [3] we report the plot of the 
parametric curve (C/(/3), {pshi{p))p) obtained by varying /3 both in the positive and in the negative region of the real 
axis. 

This curve is then compared with measures of {psui{p))E computed from molecular dynamics simulations in the 
microcanonical ensemble at different values of the energy E (Fig. [3]). Such a comparison clearly shows that the results 
obtained in the two different ensembles are identical, a transparent evidence that the equivalence of ensemble already 
exists for this system quite far from the thermodynamic limit [N = 100). 

Fig. |3] also shows that the equipartition formula cannot be used to measure the Gibbs temperature: indeed, as 
already pointed out in Section Hill the equipartition theorem can fail if the density of states ui[E) vanishes. This is the 
case of our system (Fig. [Ij, where Tq = E{E)/uj{E) should diverge for E —>■ 2N{1 + e): on the other hand the results 
obtained in the canonical and in the microcanonical ensemble clearly indicate that {psm{p))E ^ 0 as A ^ 2A^(1 + e). 


C. Spatial coherence 


In analogy with systems of point vortices discussed in Sec. IIIEl the rotators model in Eg. (1351) possesses a spatially 
ordered phase at large values of E: this can be easily understood by noting that the density of states uj{E) vanishes 
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FIG. 3: Black line: {psm{p))i 3 vs U{P) in the canonical ensemble (Eqs. (I45II and (I46II 1 as parametric functions of /3 € (— 00 , 00 ). 
Red squares: time averages of the equipartition function in molecular dynamics simulations at fixed energy E (microcanonical 
ensemble). The values for the parameters of the model are N — 100 and e = 0.5. 


in E = Em, i.e. that there is a small number of microscopic configurations corresponding to large values of E. In 
particular, the maximum of the energy Em = 2N(1 + e) is attained by the unique microscopic state where, for every 
n, Pn = TT and — (j)n-i = that is, where all the rotators are fixed (</> = sinTr = 0) and the distance among two 
consecutive rotators is A(f> = tt. As a consequence, since 4>q = Oj a-U the particles with even index (n = 0, 2,4 ...) must 
be at 0 = 0 and the others (n = 1,3,...) in (^ = tt. At smaller values of if < Em, see Fig. |1]B, such considerations 
can be extended, yielding a very similar situation; even and odd rotators must be close, respectively, to (^ = 0 or 

(f) = TT. 

Let us note that an ordered phase exists whenever, at a given energy E, the number of corresponding configurations is 
small, i.e. when u}{E) vanishes: for instance, the clustering can also be observed at small energies, when the rotators 
accumulate around ^ = 0, in order to minimize the interaction energy, see Fig. 3] B. The sign of the Boltzmann 
temperature plays a crucial role in this context, defining the features of the coherent phase. Indeed, in analogy with 
the single-particle-momentum distribution, it is easy to show that 

p{4>i - oc exp {-/3b [1 - cos{(t)i - (jri-i)] } . (47) 

When E —>• Em or if —>• 0, the inverse temperature /3 b diverges and, depending on the sign of Pb, the distribution 
Eq. (H71) peaks around (f)i — 4>i-i = 0 oy (j)i — = tt , see Fig. |3]A. 

Let us stress that not every state with negative temperature is spatially ordered: the necessary condition is a 
small corresponding phase space volume, which implies a very high energy or, equivalently, a very small negative 
temperature. The same argument applies to small positive temperatures. Of course, if negative temperatures appear, 
they signal a reduction of phase space with increasing energy, and therefore announce a more ordered structure at 
higher energy. 


V. CONCLUSIONS 

In this paper we have given a series of arguments to support the thesis of the Boltzmann temperature Tb as a useful 
parameter to describe the statistical features of a system with many particles and short-range interactions, even when 
it takes negative values. Let us draw our conclusions with a series of remarks on the role of the negative temperature 
and some comments on recent papers. 

We have shown that the temperature Tb is the proper quantity which describes the distribution of the energy 
fluctuations in the canonical ensemble. It also enters in an immediate generalization of the Mawell-Boltzmann dis¬ 
tribution to the case of “kinetic energy” which is not a quadratic function of momentum. For a particular model we 
have also demonstrated that at small ITbI (for both positive and negative values) some kind of spatial order induced 
by interactions appears, whose qualitative traits depend upon the temperature’s sign. 

If the microcanonical entropy S{e) is a convex function, independently of the sign of Tb, there is no ambiguity in 
determining the flux of energy: it always goes from the hotter system, i.e. with smaller / 3 b to the colder one (with 
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FIG. 4: A: Probability distribution function of angular distance between two consecutive rotators at high energy E — 298.96. 
B: Probability distribution function of rotators’ positions (f) in the high energy case E = 298.96 (blue squares) and in the low 
energy case E = 6.79 (red triangles). The two maxima of the high energy distribution correspond to the clusters around 0 = 0 
and 0 = TT discussed in the text. The other parameters are N — 100 and e = 0.5. 


larger /3 _b). It should be reminded that the convexity of S{e) can be violated only for very small systems or systems 
with long range interaction, both cases being very well known examples that can violate thermodynamic requirements. 

From a physical point of view it is possible to obtain the canonical ensemble from the microcanical one only for 
large systems with short range interactions. In such a class of systems, if IV 2> 1, the S'(e) is convex and it is easy 
to obtain the equivalence of the ensembles. Such a property is a fundamental requirement to obtain equilibrium 
thermodynamics, where there is no difference between thermostatted and isolated macroscopic systems. It is worth 
emphasizing that the equivalence of the ensembles only holds if one adopts the Boltzmann definition of entropy: for 
this reason, in systems exhibiting negative temperatures, where Sb and Sq are no longer equivalent in the large N 
limit, thermodynamic can be recovered for TV —> oo only through the Boltzmann formalism. 

In systems with few components and/or with long range interactions, one can still define a canonical ensemble at a 
formal level (i.e. assume that the phase space distribution is oc and then wonder about the equivalence of the 

ensembles. However such a formal mathematical approach, in our opinion, has no physical meaning. Since in presence 
of long range interactions (or equivalently a system with N = 0(1)) it is not possible to make a clear distinction 
between the system and the reservoir, it is not possible to construct systems following a canonical distribution. For 
the same reason the question of the flux of energy among two systems appears to be meaningless in those cases. 

Following Hugh [32j |. Tb can be computed via a molecular dynamics simulation, and (at least in principle) from 
the data of an experiment. The microcanonical formula (1301) . which, in most cases, allows for a practical definition 
of Tg, can fail in systems with negative Tg, therefore, as far as we know, at variance with Tb, there is not a general 
method to compute Tq in an experiment. 

We underline that the counterexamples used in Q to support the claimed inconsistency of the use of Tb are based 
on systems with very few degrees of freedom and non convex S{e). Let us note that the system in eq. (71) of @ is 
nothing but the system considered in our Section IlYl Eq. dSSl), with = 1 and e = 0: the claimed strange behavior 
of Tb is present only if fV = 0(1). On the contrary for TV ^ 1 as a consequence of the convexity of S{e) one has 
a quite natural scenario, as discussed above. In a similar way we have shown that the consistency of Tq with the 
microcanonical formula fails for large N. 

In the microcanonical ensemble the temperature Tb is a function of the total energy E. In the canonical ensemble 
the temperature Tb is a mere property of the reservoir and does not depend on the microscopic configuration of the 
system. In Q, see Sect. 3.D, the wrong concept of temperature (in non-isolated (sub)-systems) depending upon the 
energy of the microscopic configuration, see their Eq. (31), is used to claim the inconsistency of Tb. Such confusion 
seems to be persistent, see [2^ for a discussion of the topic of the (non existing) fluctuations of temperature. 

In conclusion our analysis, that applies to a large class of systems with many degrees of freedom and short-ranged 
interactions, shows that the Boltzmann temperature has the following properties: i) it is the proper quantity ruling the 
fluctuations of energy of a sub-system; ii) it can be measured by means of time-averages of a suitable observable; iii) it 
rules the direction of the fluxes of energies between two coupled systems at different initial temperatures. About the 
Gibbs temperature, we can mention that: i) the Gibbs entropy is an adiabatic invariant (although a mathematically 
rigorous proof exists only for one-dimensional systems); ii) the microcanonical formula for equipartition in general is 
not valid therefore - at variance with Tb - a simple way to measure Tq is not available. We note that the differences 
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between Tb and Tq can survive for large N, even when the ensembles are equivalent in the thermodynamic limit. 
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